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In this work, we perforin first principles DFT calculations to investigate the interplay between 
magnetic and structural properties in Ni2MnGa. We demonstrate that the relative stability of 
austenite (cubic) and non-modulated martensite (tetragonal) phases depends critically on the mag- 
netic interactions between Mn atoms. While standard approximate DFT functionals stabilize the 
latter phase, a more accurate treatment of electronic localization and magnetism, obtained with 
DFT-I-U, suppresses the non-modulated tetragonal structure for the stoichiometric compound, in 
better agreement with the experiments. We show that the Anderson impurity model, with Mn 
atoms treated as magnetic impurities, can explain this observation and that the fine balance be- 
tween super-exchange RKKY type interactions mediated by Ni d and Ga p orbitals determines 
the equilibrium structure of the crystal. The Anderson model is also demonstrated to capture 
the effect of the number of valence electrons per unit cell on the structural properties, often used 
as an empirical parameter to tune the behavior of Ni2MnGa based alloys. Finally, we show that 
off-stoichiometric compositions with excess Mn promote transitions to a non-modulated tetragonal 
structure, in agreement with experiments. 
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I. INTRODUCTION 

Magnetic Shape Memory alloys (MSMAs) are be- 
ing widely explored for many technological applications. 
Typical representatives of this class of materials are in- 
termetallic Heusler alloys with X2YZ composition (X 
and Y being transition metals). The high temperature 
phase (Austenite) of these systems is usually cubic and 
consists of four interpenetrating fee lattices. At lower 
temperatures they often assume a tetragonal structure 
(Martensite). The transition between these phases is of 
martensitic type and is usually characterized by high re- 
versibility and tunability — "— . Being magnetic, these 
materials can show a strong coupling between struc- 
tural and magnetic transitions^TiSiiii which is very ap- 
pealing for a wealth of different technological applica- 
tions. In fact, Heusler alloys have been studied quite 
extensively for their magneto-calorioi^"— and magneto- 
mechanical2r22iiiS properties and more recently for their 
applications in energy conversions^. 

Ni2MnGa is a prototype representative of MSMAs. 
The stoichiometric alloy is ferromagnetic below a Curie 
temperature of Tc — 365 K— . The martensitic transition 
to a tetragonal phase, that takes place at Tm — 200 K, 
preserves the ferromagnetic character, albeit with a vari- 
ation in the magnetic anisotropy. The martensitic phase 
is observed to have a modulated structure, where the 
parallel (110) planes are shifted from their equilibrium 
position along the [110] direction with a period of five 
unit cells (5M) and a tetragonal distortion of c/a ~ 0.94 
^~i^. Under applied stress, an orthorhombic structure 



with a modulation of seven unit cells (7M) has also been 
observed 20,21,24,25 ^ Instead, a non-modulated, tetrago- 
nal martensite with c/a ~ 1.2 can be obtained for off- 
stoichiometric compositions at lower temperaturea^li^ 
(T < 193 K). Several groups have also claimed the 
existence of a prc-martensitic phase below a transition 
temperature of TpM ~ 260 K, where the structure is 
modulated with a period of three unit cells (3M)^. 
Co-existence of several modulated and non-modulated 
phases in a twinned structure has also been observed ex- 
perimentallj^. 

Ni2MnGa and related Heusler alloys have been stud- 
ied computationally in several works that investigated 
their electronic and structural properties. These studies, 
mostly based on the local density (LDA) or the gener- 
alized gradient approximation (GGA) for the exchange- 
correlation functional, have been successful in describ- 
ing the electronic, structural^i"— , vibrational proper- 
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, and computing the phase diagrant22i2S and the 
Curie temperature^^— for these types of compounds. 
A particularly well studied aspect is the relative stabil- 
ity of the cubic austenite and a non-modulated marten- 
site phases with c/a ^ 1, which has been investigated 
by first-principles calculations using various functionals 
^-33^ With few exceptions, these studies have yielded 
qualitatively similar results that, while in disagreement 
on the amount of the tetragonal distortion, generally 
predict a stable non modulated martensite phase with 
c/a > 1. Several mechanisms have been proposed as re- 
sponsible for the structural distortions of this material, 
such as the Jahn- Teller effecii^iiMdi and Fermi surface 
nesting, which was found to be related to phonon mode 
softening^iil^ii^. It was also found that the number of 
valence electrons per atom (e/a) can control the rela- 
tive stability of different phases, the martensitic transi- 
tion temperature and the softening of the phonon modes. 



23,47,48 _ However, the microscopic mechanisms responsi- 
ble for these effects have not yet been completely un- 
derstood. As a consequence, engineering these materials 
for specific applications relies on the empirical optimiza- 
tion of alloy composition and is mostly based on a trial 
and error approach. Experiments have also found that 
the samples in which a non-modulated martcnsitic phase 
is stable are off-stoichiomctric ones, with an excess Mn 
content in their compositions 12,23,26,49 ^ j^^ fact, a related 
compound Ni2MnSn has been shown to have a low tem- 
perature tetragonal structure with cj a > 1 only with a 
Mn content larger than the stoichiometric composition, 

both experimentally and by first-principles calculations 
35 

In this work, we perform a detailed study of the rela- 
tive stability of the cubic and non-modulated martensitic 
structures and try to clarify some of the conflicting re- 
sults in the literature. After computing the effective ex- 
change interaction parameters between Mn atoms with 
constrained density functional theory (DFT), we deter- 
mine an approximate magnetic energy using the Heisen- 
bcrg model. We show that the difference in energy be- 
tween the austenite and martensite phases can be ac- 
counted for by the variation of Heisenberg magnetic cou- 
pling. By using the DFT-I-U functional^"—, we also 
demonstrate that electronic correlations play an impor- 
tant role in determining the effective exchange interac- 
tions. In particular, DFT-I-U suppresses the stability of 
the non-modulated martensitic structure for stoichiomet- 
ric compositions. These results are interpreted using the 
Anderson impurity model^, where Mn atoms are treated 
as a periodic array of magnetic impurities embedded in 
a conduction electron gas composed of electrons on Ni d 
and Gap orbitals, which mediate super-exchange interac- 
tions between Mn atoms. Using this idea, we are able to 
predict the strength of the super-exchange interactions, 
and to show that the balance between the super-exchange 
through Ga p and Ni d orbitals determines the stability 
of the structure. We also study the effect of the number 
of electrons on the effective exchange interactions, by in- 
jecting electrons into the unit cell. We show that the 
number of injected electrons affects the Fermi surface, 
which in turn modifies the strength of super-exchange 
interactions and the stability of the martensitic phase. 
Finally, we show that the martensitic phase can be sta- 
bilized when the Mn content in the material exceeds the 
stoichiometric composition value, which is in agreement 
with the experimental observations. 

This paper is organized as follows: In section |lTl 
we summarize the computational methods we have em- 
ployed. In section Hin wc discuss calculations performed 
in the austenite and martensite phases, and discuss their 
relative stability. In section HVl we present the calcula- 
tion of effective exchange parameters between Mn atoms. 
In section|Vl we discuss the underlying magnetic coupling 
mechanisms between Mn atoms. In section IVIl we dis- 
cuss the effects of the number of electrons in the unit cell 
and the extra Mn impurities on the structural proper- 



ties. Finally in section FVlIl we propose some conclusive 
remarks. 



II. COMPUTATIONAL METHOD 

Calculations in this paper were performed using 
the plane-waves pseudopotential implementation of 
DFT contained in the "pwscf" code of the Quantum 
ESPRESSO package^. The generalized gradient ap- 
proximation (GGA) functional with the Perdew-Burke- 
Ernzherof (PBE) parametrization^ was employed for the 
exchange-correlation energy. The Ni, Mn and Ga atoms 
were all represented by ultrasoft pseudopotcntials^. The 
electronic wavefunctions and charge density were ex- 
panded up to kinetic energy cut-offs of 45 Ry and 480 Ry 
respectively. The Brillouin zone integrations were per- 
formed using a 12x 12x 12 Monkhorst and Pack special 
point grids^ and a Mcthfcssel and Paxton smearing of 
the Fermi-Dirac distribution function^, with a smearing 
width of 0.01 Ry. 

In order to improve the description of electronic 
correlation and localization, the Hubbard-model-based 
DFT-I-U corrective functional was employed. The on- 
site Coulomb repulsion parameters were computed using 
the linear response method introduced in Ref. ISOl The s 
states of Ni, Mn, and Ga were treated as a charge "reser- 
voir" as described in Ref. [6II The linear response calcu- 
lation of U provides the intersitc interaction parameters 
(y)^ as well; however, we found that they are signif- 
icantly smaller than the on-site interaction parameters 
(L'^). therefore we neglect them in the calculations pre- 
sented in this paper. The crystal structures presented 
in this paper were generated using the XCrysDen soft- 



war( 



,63 



III. GGA AND GGA+U CALCULATIONS 

A. Cubic phase 

The high-temperature {T > Tm) cubic phase of 
Ni2MnGa has L2i structure where the Ni, Mn and Ga 
atoms form four inter-penetrating fee lattices (two of 
which are occupied by Ni), as shown in Fig. [TJ Also 
shown in Fig. [TJ is the tetragonal unit cell, whose basal 
unit vectors are parallel to [110] directions of the L2i 
structure. To describe the tetragonal deformation of the 
unit cell, the ratio c/a is generally used as measured from 
the cubic cell (c/a = 1). The equilibrium lattice param- 
eter and magnetization determined with GGA calcula- 
tions reported in Table |T] are in good agreement with 
experiments^. Our calculations show that the Mn d 
states split into two groups of three and two degener- 
ate states, that according to tetrahedral group notation, 
arc indicated by t2 and e respectively, whose density of 
states, obtained with the GGA functional are shown in 
Fig. [21 This shows that the Mn d states are subject 
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FIG. 1: (Color online) The L2\ structure is shown on the left, 
and the tetragonal unit cell is shown on the right. Both unit 
cells correspond to the cubic phase. 
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FIG. 2: (Color online) The projected density of states calcu- 
lated using the GGA functional. The Fermi level is set to the 
zero of energy. The inset shows the density of states in the 
close proximity of the Fermi level. 



to a cubic crystal field. More precisely, with respect 
to the Cartesian coordinates on the tetragonal unit cell 
(when cja = 1) in Fig. [1] dxy and d^^i orbitals of Mn 
are the e states and d^xi dzy, dx^_y2 orbitals are the ^2 
states. Similarly, the Ni d orbitals arc also split into e 
and ^2 states, with dxy and dz^ orbitals being the e and 
dzx, dzy, dx2-y2 orbitals being the ^2 states. 

Another feature clearly visible in the density of states is 
that the magnetization density is mainly localized on the 
Mn d orbitals. The peaks corresponding to majority and 
minority states are clearly separated, with minority spin 
states nearly empty. The localized nature of magnetiza- 



TABLE I: Calculated lattice parameter and magnetizations 
of each atom in the cubic phase. 

ao (A) flMni^s) /^Ni (ms) MGa(/is) /itot (/is/cell) 



GGA 5.83 
GGA+U 5.83" 



3.67 
4.52 



0.34 
0.16 



-0.13 
-0.13 



4.22 
4.80 



tion on the Mn d orbitals suggest that on-site Coulomb 
repulsion could play a dominant role. It can also be ex- 
pected that GGA functionals would not be able to ac- 
count for this interaction accurately, since they tend to 
delocalize the electronic density. To improve the descrip- 
tion of localized electrons responsible for magnetism, we 
adopted the DFT-hU approach ^-^. 

The GGA+U calculations were performed by using 
U = 5.97 eV on the d states of Mn. The value of the U 
correction is calculated from the linear response approach 
of Ref [50, in a 16 atom super cell (four times larger than 
the unit cell). We have explicitly verified that the cal- 
culated Hubbard parameters are unchanged when larger 
supercells are used, as explained in Ref. [50. Although 
U = 5.97 eV might seem unexpectedly large for a metal- 
lic system, it is used for describing localized electrons on 
Mn d states, which have vanishing contribution to the 
density of states at the Fermi level. On the contrary, the 
metallic character is mainly due to Ni d (e) and Ga p 
states that dominate the density of states at the Fermi 
level. The t/ on Ni d and Ga p states, as well as the in- 
teraction terms between electrons in different manifolds 
(e.g., d and s states of Mn) or on different sites (inter- 
site interactions V—) are neglected in this work as they 
are found to be largely irrelevant for the conclusions of 
this work. Also, the calculation of the electronic parame- 
ters was not performed self-consistently as in Ref. l65l nor 
they were recomputed for different distortions of the unit 
cell^, as their variation with the crystal structure was 
found to be unimportant. The density of states obtained 
from the GGA-(-U functional is shown in Fig. [31 The most 
significant difference from the density of states calculated 
with the GGA functional in Fig. [5J is the increase in the 
splitting between the majority and minority spin peaks of 
the Mn d states (consistent with other DFT+U calcula- 
tions), and their complete disappearance from the Fermi 
level. At this point, we would also like to comment on the 
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FIG. 3: (Color online) The projected density of states calcu- 
lated using the GGA-I-U functional. The Fermi level is set to 
the zero of the energy. The inset shows the density of states 
in the close proximity of the Fermi level. 

photo-emission type experiments performed in the litera- 



ture, which could be compared to the calculated density 
of states. These studies have become available relatively 
recently on Ni2MnGa(Sn)^i^'^ and the data presented 
were able to provide information mainly on states close 
to the Fermi level. These states are mainly of Ni e type, 
and as Refs. I35ll45l also suggest, they are responsible for 
the nesting of the Fermi surface. Our GGA+U calcula- 
tions modify almost entirely Mn d states, pushing them 
to lower energies, where experiments are not conclusive. 
Instead, the U correction has little effect on the states at 
the Fermi level, and therefore our results are still consis- 
tent with existing data. 

The more complete filling of electrons on the major- 
ity spin states of Mn leads to a larger magnetization, 
4.80 ^B per unit cell. The magnetization of each atom, 
computed with the GGA-I-U functional, is shown in Ta- 
ble m In our calculations, we fixed the cubic lattice pa- 
rameter to the value obtained from the GGA functional 
for simplicity. As can be seen from Table HI while the 
magnetization of Ga remains at the same value, the mag- 
netization of Mn increases and the magnetization of Ni 
decreases compared to GGA. The GGA+U results are 
in overall agreement with the neutron scattering data^, 
even if the Mn magnetization is slightly overestimated. 
In this regards, it should be noted that the extraction 
of the magnetic moment of each atom from the data re- 
quires a priori assumptions, which affect the overall out- 
come of the analysis, and leads to large uncertainties. 
We can obtain a better understanding of the effects of 
the Hubbard correction, by considering the trace of the 
occupation matrices n^^^, defined in Ref. |50|, which is a 
measure of how many electrons localize on a site / with 
spin a (while n^ = ^Ylim a '^Inm. i^ '^^ total number of 
electrons on site /). In the GGA ground state, we found 
that nj^i-d — 8.83, nun-d — 5.47 and noa-p — 2.43 
while in GGA-I-U, UNi-d — 8.90, nun-d — 5.21 and 
UQa-p — 2.46. In the GGA ground state, there are 



^Mn-d 



0.9 minority spin electrons shared between the 

i2 states of Mn, compared to n\f^_^ ~ 4.57 of major- 
ity spin. In GGA+U, which favors integer occupation 
of orbitals, the minority spin d states are almost empty 



with n\j^^ 



0.34, while "n-Mn-d — 4.86 for majority 



spin, closer to being fully occupied. While some of the 
minority spin electrons in Mn ^2 states are transferred 
to the majority spin states in GGA+U, the rest mostly 
contribute to the conduction electrons of Ni d type (since 
increase in Ga p occupation in GGA+U is small). This 
electronic re-organization in the GGA+U ground state 
has important consequences in the structural properties 
of the material, as shown in the next section. 



B. Tetragonal distortions 

We first discuss the electronic structure of the low- 
temperature phase of Ni2MnGa, with an experimental 
tetragonal distortion of c/a ~ 1. 2^^'^^ . This phase has 



often been studied in the literature for the stoichio- 
metric compound, as many calculations obtain it as a 
ground state. In order to investigate the electronic re- 
organizations when the unit cell is tetragonally distorted 
at constant volume, we have computed the density of 
states using both the GGA (shown in Fig. S]) and the 
GGA+U (shown in Fig.© at c/a = 1.2. The value of U 




FIG. 4: (Color online) The projected density of states calcu- 
lated using the GGA functional at c/a — 1.2. The Fermi level 
is set to the zero of energy. The inset shows the density of 
states in the close proximity of the Fermi level. 
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FIG. 5: (Color online) The projected density of states calcu- 
lated using the GGA+U functional at c/a = 1.2. The Fermi 
level is set to the zero of energy. The inset shows the density 
of states in the close proximity of the Fermi level. 

was set to the calculated value from the cubic phase. We 
have verified the validity of this approximation by recom- 
puting U for various tetragonal distortions (at constant 
volume), and we found that U is almost constant. 

In the tetragonal cell, Mn and Ni t2 and e orbitals split 
(the ^2 orbitals split into a doublet d^x, dzy and a sin- 
glet dx2-y2, while the e splits into two singlets) as can 
be seen from Figs. [5] and [51 Some previous works in the 
literature relate the existence of the tetragonal phase at 
c/a ~ 1.2 to a Jahn- Teller type distortion^i^i^. From 



the inset in Fig. |4l we observe that the Ni minority spin e 
peak, which represents the dominant contribution at the 
Fermi level splits, forming a shoulder towards the lower 
energy direction. This shoulder is due to the lower en- 
ergy (ij,2 orbitals, while dxy orbitals are pushed to higher 
energy. However, peaks corresponding to both types 
of orbitals move above the Fermi level while maintain- 
ing a tail below it. This feature is much stronger in 
GGA-fU, where both minority spin Ni e orbitals move 
entirely above the Fermi level, as shown in the inset of 
Fig. [5] Jahn- Teller type distortions would tend to place 
the Fermi level between the split orbitals, resulting in 
a fully occupied lower energy and an empty higher en- 
ergy orbital. Instead, we observe that the split orbitals 
both move above the Fermi level. Therefore, we argue 
that another mechanism is responsible for determining 
the energy landscape of the material and in particular, 
the stability of the tetragonal phase. In order to investi- 
gate this point further, we have performed total energy 
calculations in GGA and GGA-I-U as a function of the 
tetragonal distortion {c/ a) at constant volume, which are 
shown in Fig. [HI As can be seen, the GGA functional 
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compound Ni2Mni.25Gao. 75. Instead, Refs . I28l l 30ll33 have 
found that an energy minimum for c/a < 1 is only pos- 
sible for modulated structures. Experimentally, a distor- 
tion with c/a < 1 seems to correlate with the modulation 
of the structure 20,22,24,25 ^ j^ tetragonal non-modulated 
phase seems to be stable only for significant deviations 
from the stoichiometric compositionsi^i^i^^iiS. There- 
fore, experiments seem to not support the existence of 
a non-modulated tetragonal phase at the stoichiometric 
composition, which is in contrast with the result of the 
GGA calculation. 

In GGA-I-U, the global energy minimum at c/a ~ 1.23 
is suppressed as shown in Fig. [S] This finding indicates 
that obtaining a more accurate description of electronic 
localization on d states is important to predict the rel- 
ative stability of different phases. The suppression of 
the stable tetragonal phase shows that GGA+U is in 
better agreement with experiments, compared to GGA. 
In order to understand the relative energies of different 
phases, and their relation to magnetism in GGA and 
GGA-I-U, we have computed the magnetizations of each 
atom as a function of c/a, and the results arc reported in 
Figs. [71 [51 [ni An interesting feature that is apparent from 




FIG. 6: (Color online) Energy of the ground state per unit cell 
as a function of c/a at constant volume in GGA and GGA-f U. 
The zero of energy is set to be at c/a = 1 in each case. 



FIG. 7: (Color online) Magnetization of Ni atoms as a func- 
tion of c/a computed using GGA and GGA+U. 



yields a local minimum at c/a = 1 and a global mini- 
mum at c/a ~ 1.23. The energy difference between the 
two minima is very small and is around 4meV. Numer- 
ous studies in the literature have previously found similar 
results to our GGA calculation—"—!^. However, these 
studies have identified energy minima at slightly differ- 
ent values of c/a, and proposed different explanations 
about the relative stability of different phases. For ex- 
ample, Ref. [29| has identified a local energy minimum for 
a non- modulated orthorhombic structure at c/a ~ 1.11 
and b/a ~ 1.04, besides the local minimum at c/a ^ 1 
and a global minimum at c/a ~ 1.2. Ref. [23 has found 
a local energy minimum at c/a ~ 0.94 and a global 
energy minimum at c/a ~ 1.3 for the non-modulated 
structure, while Ref. [33| found a local energy minimum 
at c/a ~ 0.94 for the off-stoichiometric non-modulated 



Fig. [7] is that the energy minima in the GGA functional 
corresponds approximately to extrema of magnetization 
of the Ni atom (i.e. the minimum magnetization of Ni is 
obtained at c/a = 1, while the maximum magnetization 
is obtained at c/a ~ 1.2). On the other hand, the mag- 
nitude of the magnetization of Ga atoms has the lowest 
value around 1.1 ^ c/a < 1.2 but increases monotonically 
for c/a > 1.2. The Mn atoms have the highest magne- 
tization at c/a ~ 1, and it decreases monotonically for 
c/a > 1.2, with a rather featureless profile. Therefore, 
the magnetizations of Ni and Ga atoms seem to corre- 
late with opposite minima in the total energy as a func- 
tion of c/a. The oscillatory variation of magnetization 
of Ni atoms suggests a correlation with electronic screen- 
ing, since the dominant contribution to the Fermi surface 
comes from the minority spin electrons of Ni atoms. It 
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FIG. 8: (Color online) Magnetization of Mn atoms as a func- 
tion of c/a computed using GGA and GGA+U. The insets 
show a close up view of the magnetizations separately. 




FIG. 9: (Color online) Magnetization of Ga atoms as a func- 
tion of c/a computed using GGA and GGA-I-U. 



is well known that, when an impurity is embedded in a 
gas of conduction electrons, its potential is screened and 
leads to an oscillatory variation of the charge density of 
free electrons, whose period depends on the Fermi mo- 
mentum (Friedel oscillations). If Mn atoms are treated 
as a periodic array of "impurities" and the Ni d states as 
the conduction electron sea, a similar oscillatory behav- 
ior can be expected. Since the minority spin e electrons 
of Ni provide the dominant contribution to the density of 
states at the Fermi level, they contribute to the screening 
of the impurity potential most significantly. According to 
this picture, the oscillations in the magnetization of Ni d 
states can be thought of as due to spin-polarized Friedel 
oscillations. Notice that the oscillations arc with respect 
to c/a and not the distance between the Ni atoms. As 
c/a increases, the distance between Ni atoms decreases 
in the basal plane, while the vertical distance between 
them increases, which introduce a more complicated de- 
pendence of the oscillations on c/a than in the free elec- 
tron case (through the anisotropics of the Fermi surface). 



We therefore treat these oscillations at a qualitative level. 



IV. EXCHANGE INTERACTION BETWEEN Mn 
ATOMS 

The ferromagnetic ordering is a consequence of the 
magnetic interactions in the system, mainly between Mn 
atoms. Therefore, we can approximate the magnetic in- 
teraction energy through the Heisenberg Hamiltonian, 
that depends only on the magnetization of Mn atoms: 



'l-mae 






(1) 



In Eq.(IT]) Jij are the exchange parameters. Si's are mag- 
netizations of Mn atoms and (i, j) denotes that the sum 
runs over nearest neighbor Mn atoms. This approxima- 
tion is justified by our expectation that the strength of 
the exchange parameter decays rapidly with distance. In 
order to compute the exchange parameters J^.j, we use 
constrained DFT calculations in a 16 atom supercell, that 
contains four Mn atoms. In this calculation, the magnetic 
moment of one of the Mn atoms is flipped to an opposite 
direction with respect to the rest of the Mn atoms (while 
its magnitude turns out to remain unchanged). The en- 
ergy difference with the original ferromagnetic configura- 
tion is then mapped on the Heisenberg model of Eq.([l]) 
to obtain Jij. This calculation is performed for a range 




FIG. 10: (Color online) Mn atoms in the supercell used to 
compute the exchange parameters. 

of c/a values. For the tetragonal cells, there arc two rel- 
evant exchange parameters; one in the x — y plane ( Jf,) 
between the Mn atoms connected in the [110] direction 
and one in the y — z plane ( Ji,) between the Mn atoms 
connected in the [Oil] direction as shown in Fig[TUl Nat- 
urally, the two exchange parameters J^ and Jb are equal 
in the cubic limit {Jy = Jb = ^cubic at c/a = 1). The 
magnetic energy per supercell is given by 



"Hn 



(4Jfc + 8J„) {Sl + S+S-] 



(2) 



where S+ is the magnetization of the spin up Mn atoms 
and S- is the magnetization of the flipped spin Mn atom. 
The calculation is performed by first fixing the basal 
plane dimension a (thus, Jb = Jcubic) and varying c to de- 
termine J„ as a function of c using Eq. ([2]) . Then, c is fixed 
(thus, Jv = Jcubic) and a is varied to determine Jb as a 
function of a. The values of c and a for each constrained 
DFT calculation are chosen such that a? c = constant and 
therefore both J„ and Jb can be expressed as a function 
of c/a (corresponding to constant volume deformations) 
only. The resulting values of J„ and Jb as a function of 
c/a are shown in Fig. [TT] for both the GGA and GGA+U 
functionals. As can be seen from Fig. Ill[ the magni- 
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FIG. 11: (Color online) Exchange parameters J„ and Jj, com- 
puted using GGA and GGA+U functionals. 

tude of the exchange parameter Jb increases with c/a, 
while that of Jv decreases until c/a ^ 1.1 after which it 
plateaus. This is expected, since for larger c/a the dis- 
tance between Mn atoms in the x — y plane decreases 
resulting in a stronger magnetic interaction in the plane 
( Jf,), while the distance between Mn atoms in the y — z 
and x—z plane increases, resulting in a weaker interaction 
( Ju). The plateau in J„ is probably related to the verti- 
cal alignment of Mn atoms that partially compensates for 
the increased distance between (001) planes. The depen- 
dence shown in Fig. [TT] also justifies the energy minimum 
at c/a > 1 for GGA. The stronger increase in magnitude 
of Jb compared to Jy indicates that the magnitude of the 
magnetic energy in Eq.([2]) increases with c/a (since the 
increase in absolute magnetic energy due to Jb dominates 
the decrease due to J^,). Instead, for c/a < 1, the increase 
in Jb is almost symmetrically cancelled by a decrease in 
Jy so no energy minimum is observed. 

Another important result we obtain is that the mag- 
nitude of both interaction parameters are significantly 
smaller in GGA+U than in GGA. In the simple Hub- 
bard model, the exchange interaction between electrons 
localized on neighboring atomic sites can be determined 
from 2nd order perturbation theory as J ~ t^ /U , where 
t is the hopping amplitude between sites^i^. In light 
of this fact, smaller exchange parameters between Mn 



atoms in GGA+U are expected, since the effective on- 
site Goulomb repulsion is larger. Notice that the J in 
the simple Hubbard model is positive indicating an anti- 
ferromagnetic interaction, while in our constrained DFT 
calculations we find negative J's, resulting in a ferromag- 
netic ground state. This discrepancy is due to a different 
type of mechanism that leads to magnetism in our sys- 
tem compared to the simple Hubbard model. We will dis- 
cuss this point in detail in the next section. Notice also 
that both exchange parameters Jy and Jb have an oscil- 
latory modulation, as was the case for the magnetization 
of Ni d states shown in Fig. [T] This modulation seems 
to suggest that Ni d states contribute to the mediation 
of magnetic interactions between Mn atoms. Indeed, the 
oscillation of the exchange parameters between localized 
magnetic impurities in a free electron gas is a well known 
effecli^. We postpone the study of possible mechanisms 
until the next section, where we will provide a detailed 
discussion about the nature of the magnetic interaction 
between the Mn atoms. In the rest of this section, we 
study the effect of the magnetic energy on the structural 
properties of the material. Since we have computed 
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FIG. 12: (Color online) Comparison of the total ground state 
energy and the magnetic energy per unit cell in GGA. 

the parameters J„ and Jh, the magnetic Heisenberg en- 
ergy in Eq.(l2|) can be determined as a function of c/a 
(in the fully ferromagnetic case with S+ = S-). Doing 
so, we have a rough estimate of the magnetic energy in 
the system. To understand its contribution to the total 
energy, we subtract Hmag from the total energy at each 
c/a value. The result is shown in Figs. [T^] and [TJ] for 
GGA and GGA+U respectively. As can be seen from 
Fig. 1121 when the magnetic energy is subtracted from 
the total energy, the minimum at c/a ~ 1.2 disappears 
in GGA. In GGA+U, the subtracted magnetic energy 
does not change the overall profile of the dependence on 
c/a. Notice that the magnetic energy is larger in mag- 
nitude in GGA than in GGA+U. Although the magne- 
tization of Mn atoms is higher in GGA+U, the effective 
exchange parameters are significantly smaller in magni- 
tude (as seen in Fig. [TTj) , leading to a smaller absolute 
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FIG. 13: (Color online) Comparison of the total ground state 
energy and the magnetic energy per unit cell in GGA+U. 



magnetic energy. There is always an upper bound in Mn 
magnetization, whieh cannot be exceeded no matter how 
large Hubbard U is. However, larger Hubbard U will 
always lower |J| (in consistence with the perturbation 
expansion J ~ t'^/U), thus leading to a lower absolute 
magnetic energy. Therefore, we argue that the appear- 
ance of the stable structure at c/a ~ 1.2 is determined by 
magnetic energy. Then, the magnitude of the exchange 
interactions Jy and Jb needs to be determined accurately 
for a correct description of the structural properties. In 
fact, the overestimation of the exchange parameters leads 
to the (spurious) structural minimum at c/a ~ 1.2 for the 
stoichiometric compound with GGA, which is eliminated 
with the Hubbard correction. 



V. MAGNETIC COUPLING MECHANISMS 

Due to the large separation between Mn atoms in the 
crystal, it is very unlikely that a direct exchange interac- 
tion between Mn atoms could be responsible for the mag- 
netic interactions. Instead, Ni and Ga atoms between 
Mn atoms are more likely to mediate super-exchange^ 
type interactions. Indeed, hints of this possibility were 
highlighted in the previous section. For instance, we had 
found that the oscillations in the effective exchange pa- 
rameters between Mn atoms show strong similarity to 
the oscillations observed in the Ni d state magnetization. 
Magnetic interactions mediated by conduction electrons 
(like Ni d and Ga p states) have long been known, char- 
acterized by an effective exchange parameter having an 
oscillatory behavior, known as Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interactioii^a. The RKKY interaction 
arises from the polarization of free electrons in response 
to the presence of a magnetic impurity which, reaching 
other impurities, results in a magnetic interaction be- 
tween them. The oscillatory behavior in the magnetic 
interaction strength can be interpreted as result of elec- 
tronic screening. The free electron response to a pertur- 



bation is measured by the Lindhard function, which has 
an oscillatory behavior in real space and decays as the 
third power of distance (i.e. ~ cos{kp R)/{kpR)^, where 
kp is the Fermi momentum) between the impurities. If 
localized Mn d electrons are treated as a periodic array of 
"magnetic impurities" , and the Ni d and Ga p electrons 
are treated as conduction electrons, then a magnetic in- 
teraction of RKKY type can be expected. Indeed, this 
possibility has been considered in the literature before 



Before discussing the relevance of RKKY interactions 
for our system, we first provide a brief summary of the 
results from a model Hamiltonian approach, which brings 
some insight into the problem. The RKKY interac- 
tion can be derived starting from the Anderson impurity 
mode l^^i^^ , which can be expressed as 

■H = y^eic^cr +y^£' {fidt + ndi) + y^^Undtndi 

k,(T d d 

+ E (^^kcLcd.+h.c.) (3) 



k,d.i 



where ek denote the energy levels and Ukc 



cLck<T is 

the number operator for free electrons respectively (or 
Bloch electrons in the context of a periodic system). E 
denotes the energy of electrons localized on the impu- 
rity, U is the Coulomb repulsion between two electrons 
localized on the same impurity and fid t{i) ^^^ the num- 
ber operators. V^k is the hopping amplitude of electrons 
from the impurity site to the free electron states, and so 
it accounts for interactions. In writing the Hamiltonian 
(|3]), we have ignored the direct coupling between impurity 
sites: ^^ d'^^dd' cr^^Cd'a + h.c.). In fact, we assume that 
such interactions are negligible, due to spatial separation 
between Mn sites. Moreover, it was also shown in Ref. iTa 
that such interactions are generally anti-ferromagnetic 
for two impurities, when the Fermi level well separates 
the impurity energies E and E + U. This is the case for 
our system as can be seen from Figs. [2] and [3l For both 
GGA and GGA-I-U, the filled and empty d-bands of Mn 
are well separated. The simplest realization of Eq. ([3]) is 
when there is only a single impurity (d ~ 1). In this case, 
one can derive an effective exchange interaction parame- 
ter between the impurity states and conduction electrons 
in second order perturbation theory as^^ 



Jdl 



2\Vdk\^U 
\E\{U-\E\) 



(4) 



where it is assumed that the impurities are maximally 
polarized (n^^) ~ 1, (n^j,) ~ 0, £■ < 0, E+U > and the 
Fermi energy is set to zero. Similarly, the magnetization 
of free electrons, as a response to the magnetization of 
the impurity, can be derived as^ 



Mk 



1 IT. ,2 dp, 
^\Vdk\ -T-ln 



de 



E' 



{E + U)^ + A^ 



(5) 



where it is assumed that the density of free electron 
states, p(e), is a slowly varying function of e. In the 



above equation, A 



l^dkl^) p(e) is a measure of the 



mixing between free electron and impurity states (More 
precisely, A corresponds to the line width of impurity 
states, due to interactions with the free electrons). Com- 
paring equations (|4]) and ([5]) we see that the exchange 
interaction between the impurity and conduction elec- 
trons are proportional to the same factor |Vrfkp! so the 
larger the |VdkP, the larger the magnetization of free 
electrons would be. Moreover, the interaction between 
the conduction electrons and the impurity states is anti- 
ferromagnetic (Jdk > 0). The RKKY interaction can 
be obtained from a fourth order perturbation theory 
starting from a two impurity Hamiltonian analogous to 
Eq.(|3]). The resulting impurity- impurity interaction en- 
ergy is given by22 



Udd = -Jddf{kFR)Si ■ S2 



where 



Jdd ~ mkp \Jdk\ 



fix) 



2x cos(2a::) — sin(2a 



(6) 



(7) 



where m is the electron mass (replaced with effective 
mass m* in a crystal), R is the distance between im- 
purities and Si^2 are their spins. For kp R <^ 1, the 
interaction between the impurities is ferromagnetic, and 
its strength decays with the third power of their distance. 
In Ni2MnGa, we treat the Ni d and Ga p states as the 
conduction electrons. Notice that Ga p states are polar- 
ized anti-fcrromagnctically with respect to Mn d states, 
which is in agreement with the picture Anderson model 
provides. On the other hand, Ni d states are ferromagnet- 
ically ordered. This is probably related to the fact that 
the above picture does not take into account the topol- 
ogy of the Fermi surface, which is more important for Ni 
d states. Indeed, we observe from Eq.([S]) that the sign of 
the free electron polarization can change depending on 
the derivative of the density of states. 

In an alloy, the situation is more complicated than 
a model with impurities embedded in free electrons for 
mainly two reasons; first, the "impurities" i.e. localized 
Mn d states form a lattice and second, the conduction 
electrons that mediate the magnetic interactions are not 
free and therefore the strength of the interaction depends 
on the nontrivial topology of the Fermi surface. However, 
several studies in the literature have found (using a frozen 
magnon approach) the oscillatory behaviour of the Mn- 
Mn exchange parameter indicated in Eq.([7]) 40,4i^73,74 ^ 
Indeed, the possibility that the magnetic interactions are 
mediated by "X" and "Z" atoms of generic X2YZ Heusler 
alloys was considered long before in the litcrature^I, by 
studying the features in the density of states. In a simi- 
lar way, we will use some basic arguments to understand 
the nature of the magnetic interactions in our system, 
in light of the Anderson impurity model and of RKKY 
interaction. A more rigorous study would require the 
calculation of the Lindhard susceptibility using a linear 
response approachSi^S. Comparing Eqs.Q,® and (O, 
we observe that the strength of the interaction between 



the impurities is proportional to the square of the mag- 
netization of conduction electrons. Using this result, we 
model the exchange parameters Jy and Jf, in order to 
investigate the relative strength of RKKY interactions 
mediated by Ni d and Ga p states separately as 



rmodcl 



Jl 



b 
model 



{jib Mm + hb mL) ^b ^ 
{jiv Mm + hv fJ-ha) r~^ 



(8) 



where r^, and r^ are the relevant distances between Mn 
atoms (more precisely, their ratio to the cubic lattice pa- 
rameter) and the coefficients ji^2b and ji,2u arc constants 
to be determined by fitting to the data obtained from 
GGA and GGA-I-U calculations in the magnetic energy 
given in Eq.([2]). The functional form of Eqs.® does not 
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FIG. 14: (Color online) Comparison of the computed and 
fitted magnetic Heisenberg energies in GGA. 
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FIG. 15: (Color online) Comparison of the computed and 
fitted magnetic Heisenberg energies in GGA+U. 

take into account effects coming from the Fermi surface 
topology that determines the (c/a-dependent) screening 
by conduction electrons (especially from Ni d states). In- 
stead, they provide an average dependence on the mag- 
netic moments. In order to find the parameters ji,2h 



10 



and ji,2v, we fit the expressions in Eq.® using Eq.([2]) 
(^modei) to the Heisenberg energy obtained from con- 
strained DFT calculations for every c/a. The fit is ob- 
tained by the least squares method and the values ob- 
tained for ji.2b and ji.2v are reported in Table. HIl The 
comparison of the fitted and the computed magnetic en- 
ergies is shown in Figs. [T3] and [121 The quality of the 
fits can be tested quantitatively by evaluating the ra- 



tio X 



r/2 



Ei [Eme^gii) - ^modcl(«)] / E, ^mag(*): ^hcrC 



i runs over the c/a values considered. For an accu- 
rate fit, x^ < 1. We find that x^ ~ 1.1 x 10"^ for 
GGA and x^ ~ 2.5 x 10"^ for GGA-I-U. As can be 
seen from Table, [ill GGA+U predicts larger exchange 
parameters (in magnitude) for the Mn-Ni interaction 
{jiv,b) than from GGA, and smaller ones for the Mn- 
Ga interaction {J2v,b)- The larger Mn-Ga interaction 
in GGA can actually be understood using Eq.Q and 
the density of states. From Fig. [51 we observe that 
the lowest energy filled Mn d states are peaked around 
3.6 eV below Ef, while the unoccupied Mn d states are 
peaked around 1.6 et^ above Ef. This corresponds to 
\E\ ~ 3.6 eV and -\E\ + U^ff c^; 1.6 cV (where Ua« is 
the effective Coulomb interaction on Mn site) and there- 
fore from Eq.gJ we find Jun-Ga ^ 1.8 |VA,/„-Ga|^ in 
GGA. Instead, from Fig. [31 we see that the occupied 
Mn d states are peaked around — 6.4eV, while the un- 
occupied Mn d states are peaked around 4eV, which 
leads to |^| ~ QAeV and -\E\ + U^s ~ 4eV. There- 
fore, JMn-Ga ^ 0.81 \VMn-Ga\^ in GGA+U. This valuc 
is further lowered by the fact that |VA/n-Ga| in GGA-I-U 
is smaller due to stronger localization of electrons on Mn 
sites. Instead, the larger exchange parameters for Mn-Ni 
interaction in GGA-I-U is related to the factor kp f{kp R) 
in Eq.([7]). Since Ni d states have the dominant contribu- 
tion at the Fermi level, the dependence on the Fermi mo- 
mentum (to be more precise, the Fermi surface) is more 
important for the Mn-Ni exchange interaction. However 
Eq. (Hj) is not sufficient to account for the topology of the 
Fermi surface, and therefore we can not use it to estimate 
the relative strength of Mn-Ni exchange parameters, un- 
like we did for Mn-Ga exchange parameter. Instead, we 
interpret Mn-Ni exchange parameters jiv,b as describ- 
ing the magnetic interaction in an average way. The 
inaccuracy to model the Fermi surface topology in the 
magnetic interactions, also explain why the fitted mag- 
netic energies shown in Figs [TJl and [TSl arc less accurate 
for c/a <, 1.1 but more accurate for c/a ^1.1. As can 
be seen from Fig. [71 the magnetic energy is correlated 



TABLE II: Fitted parameters in the model for J„ and Jb 
found in GGA and GGA-hU. 





GGA [Ry/^l] 


GGA-HU [Ry/^|] 


jib 


-1.94 X 10"" 


-2.73 X 10"* 


J2b 


5.63 X 10"* 


-1.61 X 10"" 


jlv 


4.99 X 10"^ 


1.68 X 19"* 


J2v 


-1.19 X 10"=* 


-3.17 X 10"'' 



with the Ga magnetization for c/a >, 1.1. However for 
c/a < 1.1, Ga magnetization is smaller and Ni d states 
start to have larger influence on the exchange parameters 
and due to the dependence on the Fermi surface topology, 
the inaccuracy in determining Jmu-Ni becomes more rel- 
evant. Another trend we observe is that the exchange pa- 
rameters |jiu,6| are in general smaller than |j2u,b|. The 
only exception is \J2b\ < \jib\ in GGA-I-U. These results 
indicate that the stable structure of Ni2MnGa is deter- 
mined by a competition between the exchange parame- 
ters of Mn-Ni and Mn-Ga interactions, that determines 
the magnetic energy as a function of c/a. The model in 
Eq.® provides a simple average measure of the depen- 
dence of magnetic interactions on the magnetization of 
Ga and Ni atoms, which can be used to explain different 
mechanisms that determines the relative stability of cu- 
bic and tetragonal phases. Specifically, Eq.® provides 
an accurate estimate of the dependence of the magnetic 
energy on c/a and shows that Mn-Ni super-exchange is 
prevalent for c/a ^ 1, while Mn-Ga super-exchange is 
prevalent for c/a > 1. In Table. [Ill we report the effec- 
tive interaction strengths related to these mechanisms. 
While the Mn-Ga exchange interaction is easily realized 
through the Anderson impurity model, the Mn-Ni inter- 
action requires a careful study of the Lindhard electronic 
response function, that is closely related to the Fermi sur- 
face topology. Indeed, the softening of phonon modes, 
the dependence of the stability of the structure on the 
number of valence electrons per unit cell, and the ap- 
pearance of modulated structures have been argued to 
be connected to the Fermi surface topology 37,45,46 ^ q^j, 
study is a further confirmation of these hypothesis. 



VI. ELECTRON INJECTION AND 
OFF-STOICHIOMETRIC ALLOYS 

A. Extra electrons in the unit cell 

The correlation between the total number of valence 
electrons per atom (e/a) in the unit cell and structural 
properties of Heuslcr alloys, such as the Martensitic tran- 
sition temperature, the stability of the cubic structure 
and the vibrational spectra, has been studied extensively 
in the literature22i^iiZrilS. In order to investigate the 
effects of e/a, we have introduced fractional electrons in 
the unit cell in the range 0.1 < e < 0.7, compensating the 
extra charge with a uniform positively charged jellium. 
In each case, we have relaxed the unit cell and optimized 
the cubic cell parameter which is ao — 5.91 A for e = 0.1 
(slightly larger than the original cell) and increases mono- 
tonically up to oq ~ 6.48 A for e = 0.7. In Fig. \W\ we 
show the ground state energy (the zero of the energy cor- 
responds to the cubic phase in each case) as a function 
of c/a for a selection of extra fractional charge using the 
GGA functional. Although GGA is not sufficient to pre- 
dict the relative stability of different structural phases, it 
will be the functional of choice in this part to discrimi- 
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FIG. 16: (Color online) Energy of the unit cell as a function 
of cja for different values of added charge, using GGA. 



nate the effects of varying e/a on the exchange couplings 
from those due to the electronic localization. The Hub- 
bard (+f7) correction will be taken into account when 
computing the value of the exchange parameters in pres- 
ence of a fixed amount of extra electrons. At the level of 
the Anderson model, GGA overestimates |Vdk| m Eq.(|4]) 
and underestimates the effective on-site Coulomb repul- 
sion \J . Instead, by changing e/a we modify the Fermi 
surface directly and IJ is unaffected at first approxima- 
tion (of course, since the problem is a self consistent one, 
the Fermi surface would eventually affect it, but only at 
a higher order). Indeed, we have investigated the den- 
sity of states for ground states with extra charge, and 
found no significant difference with the original system, 
except a shift of the Fermi energy to higher values for 
larger e/a. The effect of adding charge is most dominant 
in Ni e states whose density of states is shown in Fig [171 
As can be seen, for larger e/a the peaks just below the 
Fermi level becomes stronger since more electrons need 
to be accommodated below the Fermi energy. Therefore, 
the addition of extra charge in the system modifies the 
k\< f{kF R) term in the RKKY interaction of Eq.([7]) for 
the free electron case, and an analogous averaged quan- 
tity depending on the Fermi surface corresponding, in 
our case, to model parameters jiv,b m Eq.(|S]). Indeed, as 
we show in Fig. [THl as the extra charge injected into the 
system increases, the Friedel-like oscillations (which also 
depend on the same factor kp /{kp R) for free electrons) 
have larger wavelength. Since the period of charge oscil- 
lations is inversely proportional to kp, Fig. [T51 suggests 
that the larger the charge e, the smaller kp is. Therefore, 
from Eq.(l7l), we expect a smaller exchange coupling, that 
leads to a smaller Heisenberg magnetic energy. This is 
indeed visible in the total energy calculations in Fig. [11] 
as the stable structure at c/a ~ 1.2 disappears when 
e > 0.1. In order to verify this claim, we have computed 
the exchange couplings and the magnetic energy, as we 
did in the previous section, for e = 0.5. In Fig. 1191 we 
show the Mn-Mn exchange parameters Jy and Ji, com- 



FIG. 17: (Color online) Density of states of Ni e states around 
the Fermi level. The zero of energy for each curve is set to 
the Fermi energy of the corresponding charged state, which 
increases with charge e. 
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FIG. 18: (Color online) rNid = ?lNid/?^Nid"^ a-s a function of 
c/a for a couple of extra charge. 
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FIG. 19: (Color online) Exchange parameters J„ and Jb com- 
puted using GGA and GGA-I-U functionals with additional 
e = 0.5 electrons. 



puted with constrained DFT calculations using GGA and 
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GGA+U functionals {U has the same value as in the orig- 
inal structure). Notice that, as was the case for the orig- 
inal neutral cell, the exchange parameters in GGA-I-U 
are smaller than in GGA as expected. Comparing with 
the values obtained from the original structure with no 
extra electrons, shown in Fig.[TTl we observe that the ex- 
change parameters are smaller for e = 0.5. This can also 
be seen from Figs. [501 and [5T] in comparison with Figs.fT^ 
and [131 where we show the difference between the mag- 
netic Heisenberg energy and the total ground state en- 
ergy. Moreover, the oscillatory behaviour in the original 
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FIG. 20: (Color online) Comparison of the total ground state 
energy and the magnetic energy per unit cell in GGA with 
extra charge e = 0.5. 




FIG. 21: (Color online) Comparison of the total ground state 
energy and the magnetic energy per unit cell in GGA+U with 
extra charge e = 0.5. 

e = system is also suppressed. Both these observations 
are in agreement with the previous discussion based on 
Friedel-type oscillations having higher period for e = 0.5, 
thus leading to smaller exchange parameters Jy^b with 
suppressed oscillatory behaviors. Therefore, we expect 
the model given in Eq.(IS]) for the exchange parameters 
to be more precise in case of e = 0.5, since the effects of 
Fermi surface is less pronounced. Indeed, this is what 
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FIG. 22: (Color online) Comparison of the computed and 
fitted magnetic Heisenberg energies per unit cell in GGA with 
e = 0.5. 
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FIG. 23: (Color online) Comparison of the computed and 
fitted magnetic Heisenberg energies per unit cell in GGA+U 
with e — 0.5. 



we observe in Figs. [21] and [531 In the previous section, 
we have argued that the inaccuracies between computed 
and fitted exchange parameters arise due to the model 
in Eq.((S]) being unable to capture effects of the Fermi 
surface topology. In the current case, these effects are 
suppressed, thus the model adheres more closely to the 
direct calculation. We find that x^ ~ 1.92 x 10~^ in 
GGA and x^ ^ 1-2 x IQ-^ in GGA+U, almost an order 
of magnitude smaller than the case for the neutral cell. 
In Table HID we report the values of the fitted parameters 
for Jy and Jt, which are on average smaller than the fit- 
ted values of Table [ill for the neutral cell as expected. As 
was the case for the neutral cell, the Mn-Ga interaction 
is suppressed in GGA+U compared to GGA. For Mn-Ni 
interactions, we again observe the fitting parameters to 
be an average of the oscillatory behavior, with random 
signs. Notice that the stable structure is again deter- 
mined from a competition between Mn-Ni and Mn-Ga 
super-exchange interactions, as described on average by 
Eq.dHl). 
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We would like to remark that the suppression of the 
tetragonal structure in GGA+U, as demonstrated in sec- 
tion |lll] does not preclude the stabilization of the marten- 
sitic phase for the stoichiometric compound. As we 
have shown, the relative stability of the cubic and non- 
modulated tetragonal phases for the stoichiometric com- 
pound is determined by magnetic energy. Instead, the 
appearance of modulated structures, as well as phonon 
mode softening, are expected to be related to Fermi sur- 
face nestingSi^i^. Magnetic interactions also depend 
on Fermi surface properties, as discussed in this section. 
However, the dominant contributions come from elec- 
tronic localization on Mn d states, as demonstrated in the 
previous section. Since the Hubbard U correction does 
not modify the states responsible for the Fermi surface 
nesting (Ni e states), GGA-fU is not expected to suppress 
modulated structures. Indeed, preliminary GGA-f-U cal- 
culations (not reported here) on the 5M modulated struc- 
ture revealed the existence of stable martensitic phase 
with c/a ~ 0.92, in agreement with experiments^i"— . 
Moreover, larger e/a values suppress the tetragonal phase 
as can be seen in Fig. ^E\ in analogy with the effect of 
U. At the same time it has been shown in the litera- 
ture that larger values of e/a lead to stronger softening 
of phonon modes and some of the elastic constants ^d^. 
These results, which are both based on GGA, are indica- 
tions of the fact that the two phenomena (pure tetragonal 
distortions and modulations/phonon mode softening) are 
largely independent. 



B. OfF-Stoichiometric Compounds 

The effects of extra Mn in Ga sites has been discussed 
in detail in the literature before 27,34,48,49,81 g^j-^^^ y^^^ jj^ 

place of Sn, for the related Heusler alloy Ni2MnSn in 
Ref. |3^. It has been shown that the Mn in Ga sites 
prefers an anti-ferromagnetic ordering with respect to 
normal Mn sites^L^. Moreover, the samples that un- 
dergo a martensitic transition to a tetragonal structure 
with c/a ~ 1.2 always have larger Mn content than the 
stoichiometric compounds in Ni2MnGapi^i^'2£iiS. In the 
picture we have provided in the previous sections, the 
structure with c/a ~ 1.2 is stabilized by magnetic inter- 
actions between Mn atoms. In this section, we show that 
the c/a ~ 1.2 phase appears only in Ni2Mni+j;Gai_:c 
compounds, when x ^ 0.5, using the GGA+U functional. 
Instead, the GGA functional produces inaccurate stabil- 



TABLE III: Fitted parameters in the model for J^ and Jt 
found in GGA and GGA-I-U with extra charge e = 0.5. 





GGA [Ry/^l] 


GGA-fU [Ry/^l] 


jib 


5.94 X 10"'' 


-1.64 X 10"* 


J2b 


-9.47 X 10"'' 


9.67 X lO-'^ 


jlv 


-6.66 X 10-5 


5.97 X 19-5 


J2v 


-2.67 X 10"* 


-2.30 X 10-* 



ity profile as a function of c/a for these alloys. 

The total energy vs c/a profiles for Ni2MnGa, 
Ni2Mn1.25Gao.75 and Ni2Mni.5Gao.5 using GGA and 
GGA+U are shown in Figs. [Ml and [25l For the off- 
stoichiomctric compounds, we have used a 16 atom su- 
pcrcell, and replaced one Ga with a Mn in case of 
Ni2Mn1.25Gao.75, and two Ga with Mn atoms in case 
of Ni2Mni 25Gao.75. In the latter case, the second Mn 
impurity is placed in the same x — y plane with the first 
Mn impurity. We have also studied the system where the 
second Mn impurity is placed in a different x — y plane 
with respect to the first Mn impurity and found that 
both systems have similar electronic and structural prop- 
erties. For both off-stoichiometric compounds, we have 
relaxed the unit cell, optimized the cubic lattice parame- 
ter (using GGA) and studied constant volume tetragonal 
distortions. We have used the same U on Mn atoms 
calculated for the stoichiometric cell. Although a more 
precise approach should involve a re-calculation of the 
Hubbard U in this case, as a first approximation, the 
same U value is sufficient to study the general trends. 

We have found that an antiferromagnetic ordering of 
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FIG. 24: (Color online) Energy of the ground state as a func- 
tion of c/a at constant volume per unit cell in GGA. The zero 
of energy is set to be at c/a = 1 in each case. 

extra Mn atoms with respect to the original Mn atoms 
is energetically favored in agreement with the previous 
studies^li^. This is expected, since the extra Mn im- 
purities are nearest neighbors to the original Mn atoms 
(instead of being separated by Ni and Ga atoms), and 
therefore the exchange interaction is mainly mediated by 
a direct hopping of electrons between them (i.e. through 
the term J2d,d'(^dd' cl^Cd'a + h-c), we have neglected 
previously). This type of direct hopping results in an 
anti-ferromagnetic coupling with an exchange parame- 
ter J ~ \Vdd'\'^/U. Addition of extra Mn notably in- 
creases the magnitude of magnetic Heisenberg energy, 
leading to a stronger decrease of total energy in both 
GGA and GGA+U. However, since GGA overestimates 
the magnetic couplings, it yields a local energy minimum 
at c/a < 1 and a global energy minimum at c/a > 1.2 
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FIG. 25: (Color online) Energy of the ground state as a func- 
tion of cja at constant volume per unit cell in GGA+U. The 
zero of energy is set to be at cja = 1 in each case. 



for Ni2Mn1.25Gao.75. This result is in contradiction with 
experimental data (an energy minimum for cj a < 1 was 
observed only for modulated structures). In the case of 
Ni2Mni.5Gao.5, there is no local minimum for c/a < 1.2 
and the global minimum is placed at c/a ~ 1.3 overes- 
timating the experimental value. Instead, the GGA+U 
functional yields only a global minimum for c/a ~ 1.15 
for Ni2Mni.5Gao.5. This finding is also in agreement with 
the case of Ni2Mni+a;Sni_a; where the Martensitic tran- 
sition (to a non-modulated structure) is observed only 
when X > 0.5^. 



VII. CONCLUSIONS 

In this work we have studied the electronic and struc- 
tural properties of Ni2MnGa and investigated the relative 
stability of cubic austenitc and non- modulated tetragonal 
martensite phases. We have shown that the relative sta- 
bility of these phases critically depends on the magnetic 
energy due to exchange interactions between Mn atoms. 
This result was demonstrated by modelling the magnetic 
energy with the Heisenberg Hamiltonian and computing 
the effective exchange couplings with constrained DFT 
calculations for every structure. Comparing the results 
obtained from GGA and GGA+U functionals, we showed 
that electronic localization has an important effect on the 
magnetic energy and on the relative stability of the two 
structures. In particular, GGA+U removes the spurious 
energy minimum at c/a ~ 1.2 for the stoichiometric com- 
pound. An explanation of these results was obtained by 
using the Anderson impurity model. Within this model, 
we have treated the Mn d orbitals as "magnetic impu- 
rities" embedded in conduction electrons on Ni d and 
Ga p orbitals, which mediate RKKY type interactions 



between Mn atoms. Using this theoretical framework, 
we were able to provide an explanation for the relative 
stability of cubic and non-modulated tetragonal struc- 
tures in various circumstances. We have also formulated 
a simplified model for the Mn-Mn exchange interaction 
based on a linear combination of two contributions: (1) 
super-exchange mediated by Ga p states and (2) super- 
exchange mediated by Ni d states and were able to un- 
cover the competition between these mechanisms in de- 
termining the stable structure. This model was also used 
to study the effect of e/a, and it showed that the stabil- 
ity of the non-modulated martensitic phase is suppressed 
by increasing the number of electrons. We also stud- 
ied off-stoichiometric compounds with excess Mn (sub- 
stituting Ga atoms). We have shown that, in agreement 
with experiments, the non-modulated martensitic struc- 
ture with c/a > 1 can only be stabilized by excess Mn 
content, which increases the exchange couplings J be- 
tween Mn atoms. This result was obtained using the 
GGA+U functional. The GGA functional predicts, in- 
stead, the non modulated tetragonal phase to be stable 
at any composition, in contradiction with observations. 
The study offered in this work shows that a precise ac- 
count of magnetism and magnetic interactions is essential 
to predict the relative stability of different phases and to 
rationalize the experimental observations. The simple 
RKKY interaction is able to explain most of the physics 
of the systems studied in this paper and to capture the ef- 
fects of doping. In spite of the metallic character of these 
systems, DFT+U proved to be crucial for an accurate de- 
scription of the localization of electrons and of magnetic 
properties that play a key role in the relative stability of 
different structural phases. To the best of our knowledge, 
DFT+U was not previously used to study the electronic 
structure of Ni2MnGa. However, a different Heusler al- 
loy Co2Mni_^Fe2:Si was studied with DFT+U^ and the 
importance of electronic correlations was stressed. In 
light of these results, we argue that the study of phonons, 
Fermi surface nesting and off-stoichiometric alloys should 
be revisited using the DFT+U functional. A recently de- 
veloped extension of DFPT^^ to the DFT+U scheme^i 
will greatly facilitate the development of these studies. 
Our future calculations will use the theoretical framework 
developed in this work to explore the effects of various 
substitutional impurities and to characterize modulated 
martensitic structures. 
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